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The transition form factors of B® — ► D+(2317) and B® — ► Z?+(1968) at large recoil region 
are investigated in the light cone sum rules approach, where the heavy quark effective theory 
is adopted to describe the form factors at small recoil region. With the form factors obtained, 
we carry out a detailed analysis on both the semileptonic decays B® — > £>+(1968, 2317)lDi 
and nonlcptonic decays B s — > Z?+(1968, 2317)M with M being a light meson or a charmed 
meson under the factorization approach. Our results show that the branching fraction of 
B® — > D+ (2317) fiD^ is around 2.3 x 10~ 3 , which should be detectable with ease at the 
Tcvatron and LHC. It is also found that the branching fractions of B® — > D+(l968)li>i 
are almost one order larger than those of the corresponding B® — ► Df(23l7)li>i decays. 
The consistency of predictions for B s — ► Z?+(1968, 2317)L (L denotes a light meson) in the 
factorization assumption and fcr factorization also supports the success of color transparency 
mechanism in the color allowed decay modes. Most two-charmed meson decays of B s meson 
possess quite large branching ratios that are accessible in the experiments. These channels 
are of great importance to explore the hadronic structure of charmed mesons as well as the 
nonperturbative dynamics of QCD. 

PACS numbers: 14.40.Lb, 13. 20. He, 11.55.Hx 



I. INTRODUCTION 



Enthusiasm for the open charm spectroscopy has been renewed since the announcement of 
a narrow low mass state Z? s (23±7) with unexpected and intriguing prosperities, observed in the 
D s7 r° decay mode by BaBar collaboration fl. The analysis of these charmed resonances can be 
considerably simplified in the limit of infinite heavy quark mass, when the heavy quark acts as a 
static color source so that its spin is decoupled from the total angular momentum of the residual 
light degrees of freedom. Weak production of charmed mesons in the B s meson decays induced by 
the b — > c transition serves as an ideal platform to scrutinize the KM mechanism of the standard 
model (SM), explore the dynamics of strong interactions as well as probe the signals of new physics. 
Moreover, valuable information on the inner structures of the exotic charmed mesons can also be 
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extracted from the rare decays realized via the b — ► c transition. 

On the experimental aspect, B s meson will be copiously accumulated at the LHC, which makes 
the investigations of the -B s 's static prosperities and its decay characters promising. On the theo- 
retical side, the heavy quark symmetry can put stringent constraint on the form factors responsible 
for B° -> D+ (1968, 2317) transition. As for the B s meson transitions to the lowest lying charmed 
mesons, one needs to introduce a universal Isgur-Wise function £(v ■ v'), whose normalization is 
£(v ■ v' = 1) = 1 as a consequence of the flavor conserving vector current. However, the heavy 
q uar« symmetry could not predict the normalisation of the universal form factor r l/2 responsible 
for the decays of B q meson to the doublet j£ = (0+, l+) 1/2 (2J, therefore one has to rely on some 
nonperturbative methods to deal with the B® — ► Df( 1968, 2317) transition form factors. 

Currently, there have been some studies on the semileptonic decays 5° -> Z>+(1968, 2317)/^ 



ranging from phenomenqlo gica l model [3] to QCD sum rules approach 



,y,B, 



PQCD approach 



Q] and Lattice QCD [8j, [9j, |lGj. It could be found that the available theoretical predictions vary 
from each other, hence the investigation of these modes in the framework that is well rooted in the 
quantum field theory is in demand. 

Light cone sum rule(LCSR) offers an systematic way to compute the soft contribution to the 



transition form factor almost model-independently 
dard QCD sum rule (QCDSR) technique 1(], 17, 
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131 . 1 14 , Il5| | . As a marriage of the stan- 



181 ] and the theory of hard exclusive process, 



LCSR cures the problem of QCDSR applying to the large momentum transfer by performing the 
operator product expansion (OPE) in terms of the twists of revelent operators rather than their 
dimensions Therefore, the principal discrepancy between QCDSR and LCSR consists in that 
non-perturbative vacuum condensates representing the long-distance quark and gluon interactions 
in the short-distance expansion are substituted by the light cone distribution amplitudes (LCDAs) 
describing the distribution of longitudinal momentum carried by the valence quarks of hadronic 
bound system in the expansion of transverse-distance between partons in the infinite momentum 



frame. Phenomenologically, LCSR has 
mesons and baryons in recent years 



jeen appliec 
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23|, 



widely to the investigation of the transition of 
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25, 



26. 



27|, 



28]. 



In this work, we will employ the LCSR approach to compute the B® — ► D+(1968, 2317) form 
factors, and then analyze the mentioned semileptonic modes as well as the nonleptonic decays 
B s — > D s jM, with M being a light meson or a charmed meson, under the factorization approach. 
It is expected that we can win the double benefit from such decays: gain better understanding on 
the dynamics of strong interactions and clarify the inner structures of D+(1968, 2317) mesons. 

The layout of this paper is as follows: We firstly collect the distribution amplitudes of 
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.0^(1968,2317) mesons in the section II. The equation of motion and heavy quark symmetry 
are employed to simplify the structures of hadronic wavefunctions. The sum rules for the transi- 
tion form factors B° D+( 1968, 2317) up to twist-3 are then derived in section III, where the 
relation of form factors in the heavy quark limit are found to be well respected in the LCSR 
approach. The numerical analysis of LCSR for the transition form factors at large recoil region 
are displayed in section IV. Heavy quark effective theory (HQET) is adopted to describe the 
B® — ► (1968, 2317) transitions at the small recoil region. Moreover, detailed comparisons be- 
tween the form factors obtained under various approaches are also presented here. Utilizing these 
form factors, the branching fractions of semileptonic decays B® — > D+(1968, 2317)Zi>/ and nonlep- 
tonic decays B s — > L>+(1968, 2317)M are calculated in section V. In particular, some remarks on 
the factorization of nonleptonic modes are given here. The last section is devoted to the conclusion. 



II. EFFECTIVE HAMILTONIAN AND LIGHT CONE DISTRIBUTION AMPLITUDES 
A. Effective Hamiltonian for the b quark decays 

In this subsection, we would like to collect the effective Hamiltonian for b quark decays after 
integrating out the particles including top quark, and Z bosons above scale \i = 0{m^). For 
the semileptonic b — > clvi transition, the effective Hamiltonian can be written as 

H eff (b - clvi) = ^=V cb c^(l - ts)WV(1 - 75 K (i) 

For the nonleptonic transition with AB = 1, the effective Hamiltonian is specified as 

H <® = % E X i D) ( Cl Ql+C 2 Q P 2 + Yl Ci & + C 77 Qll + Csg Q 8 g) + h.C. , 



p=u,c v «=3,...,10 



where the CKM factors are 



(2) 



\ (D) -rr T7-* J V pb V pd> for 6 ^ d transition ; 

= VpbVpD = \ W) 
I VpbV p * s , for 6 — > s transition . 

The function Qi are the local four-quark operators: 
• current-current (tree) operators 

Ql = {pb)v-A{Dp)v-A , Q P 2 = (pibj) V -A{DjPi)v-A , (4) 



4 



QCD penguin operators: 



Q 3 = (Db) v _ A ^2 q (qq)v_A, 
Q 5 = (Db) V -A^2q(qq)v+A, 



Q 4 = (Dibj) v _ A {qjQi)v-A , 

Q 6 = (Dibj)v-A {q~jqi)v+A , (5) 



electro-weak penguin operators: 



Q 7 = (Db)v-A ^q ^e q (qq) v +A , Qs = (Dibj)v-A ^q ^e q (qjqi) v +A , 

Q 9 = (Db) v _ A ^2 g ^e q (qq)v-A, Qw = (A&j)v-A ! e g(<Zj<7i)v-A , (6) 

• electromagnetic and chromomagnetic dipole operators : 

Q 7l = ^Lm b 3*^(1 + l5 )F>"b, Qs a = ^rn b Da^(l+ l5 )G^b, (7) 

where i and j are the color indices, {q\q2)v±A = <?i7^(l ± 75)92 and the sum runs over all active 
quark flavors in the effective theory, i.e., q = u,d, s, c, b. The combinations ctj of Wilson coefficients 
are defined as usual 



ai = C 2 + Ci/3, a 3 = C 3 + C 4 /3, 05 = C 5 + C 6 /3, a 7 = C 7 + C 8 /3, a 9 = C 9 + C 10 /3, 

a 2 = a + C 2 /3, a 4 = C 4 + C 3 /3, a 6 = C 6 + C B /3, a 8 = C 8 + C 7 /3, a 10 = C 10 + C 9 /3. (8) 



B. Distribution amplitudes of Z?+(1968) 



The distribution amplitudes of pseudoscalar meson Df can be defined as [3. 
{Df(P)\c(y) W s(w)\0) = -if DsP , f due^ p - k >y- ik - w ^ D (u) 



^ * ' % Jo 

(D+(P)\c(y)-t 5 s(w)\0) = -im f Ds I dwi -i<P-k)-v-*-v> ; 

J 



)(P^z u - P vZv ) \ due-^-Q-y-^Uu), (9) 



where z = y — w and it = 1 — J^p is the longitudinal momentum fraction carried by the charm 
quark. In the heavy quark limit, the chiral mass can be simplified as 



m 



m 



m Da +0(A), 



(10) 



m c + m s 

which indicates that the contribution from the distribution amplitude <j>£)(u) is suppressed by 
0(A/m£) s ) compared with that from 4> v D {u) and 4> p D {u). It can also be observed that the twist-4 
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distribution amplitude gD{u) contributes at the power of r 2 with r = j^-, therefore it can be 
safely neglected in the numerical calculations. 

In the next place, we would like to derive the relations between the distribution amplitudes 
(j> v r)(u) and (j)£)(u) in the heavy quark limit with the help of the equation of motion. Following the 



Ref. 

as 



30l | . the nonlocal matrix element with the insertion of pseudotensor current can be rewritten 



p+(P)|c( y y^ 7 5s(™)|0> = (D+(P)\c(y)w„j&8(w)\0) - i gfM „(D+(P)\c(y) l5 s(w)\0). (11) 

Differentiating both sides of the above equation with respect to w u for fj, = — and to for v = +, 
we have 



/ 



duu<p p D (u)e- i{p - k) - v ~ lkw = 0(A/m Ds ), (12) 
du[f D {u) - 4>l{u))e-^ p - k >y- ik - w = 0(A/m Da ), (13) 



with u = 1 — u. As shown in Eq. (|12p . the distribution amplitude (j? D peaks at the region of 
u ~ 0(K/m£, s ). Eq. ([13]) indicates that the distribution amplitudes 4> v D {u) and (t> p D {u) have the 
same normalizations 

^ du^Uu) = f du4> p D (u) = f du<p D (u) = Sd.. (14) 
Jo Jo Jo 

In this way, one can express the nonlocal matrix elements relevant to the pseudoscalar D s meson 

in the heavy quark limit as 

(Dt{P)\c{z)is{U)M = \J^due™ p - z ci> D {u)[ lb {JP + m Ds )] 3U (15) 

The model of 4>d(u) adopted in this work is 

M«) = /d.6«(1 - «)[1 - C D (1 - 2u)}, (16) 

where the shape parameter Co = 0.78 is determined to fit the requirement that 4>d(u) has a 

m Ds -m c 



maximum at u 



Following t 
be defined by 



C. Distribution amplitudes of _D+(2317) 



ie same philosophy, the distribution amplitudes of scalar charmed meson D* can 



<Z)+(2317)(P)|c(^5(0),|0) = \J due^i-if^miu) +m DM (I) lj $ D2 (u)} t (17) 



where D* denotes the D^(2317) meson and the normalizations of distribution amplitudes are 

•l pi 



[ du$ m (u) = f D * , [ du<S> D2 (u) = f D * (18) 
Jo Jo 

The decay constants fu* and /d* are given by 

(0\s ltl c\D* s0 (P)) = f D *P^ , (0\sc\D s0 (P)) = m D J D * Q , (19) 

where /d* = ( m c — m s)fD* s0 / r nD* s0 with m c and m s being the current masses of charm quark and 
strange quark, respectively. 

Again, with the help of equation of motion, one can find that the distribution amplitudes 
&di(u) and <&D2(u) differ at the order of A/rri£)* o ~ (jnD* ~ m c)/ m D* - Hence, for the leading 
power calculation, it is reasonable to parameterize the distribution amplitudes and ^D2{u) 

in the following form 

$Di(u) = <S>D2(u) = /d;„6u(1 -u)[1 + a(l - 2u)\ (20) 

in the heavy quark limit. /d* — (225±25)MeV has been determined from the two-point QCD sum 
rules. The shape parameter a = —0.21 is fixed under the condition that the distribution amplitudes 
^Di(u) possess the maximum at u = — -f 2 ^ — with the charm quark mass m c = 1.275GeV. It 
is worthwhile to point out that the intrinsic b dependence of the charmed meson distribution 
amplitudes has been neglected in the above analysis, which will introduce more free parameters. 

III. LIGHT CONE SUM RULES FOR FORM FACTORS 
A. Sum rules for B° — * D+(2317) transition form factors 

The hadronic matrix element involved in the B® — > D*q transition can be parameterized as 



(D*+(P)\cw 5 b\B s (P + q)) = -iif+iq 2 ^ + f^(q 2 )q,\. (21) 

sO sU 



Following the standard procedure of sum rules, the correlation function for (q ) and f D * (q ) 
is chosen as 



n„(P,<7) = ~j #xe*^D${P)\T{j^{x)Mm% ( 22 ) 

where the current j2n{x) = c( x )l^l5^( x ) describes the b — ► c weak transition and ji(0) = 
6(0)2755(0) denotes the B s channel. 
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Inserting the complete set of states between the currents in Eq. (|22p with the same quantum 
numbers as B s , we can arrive at the hadronic representation of the correlation function 

n , p x _ , (D*+(P)\c(0h^H0)\Bs(P + g))<B a (P + g)|6(0)i758(0)|0) 

. (/j; + (P)|c(0)7 M 75fc(0)l^(^ + + g)|6(0)i75«(0)|0) f23) 

^ m 2 -(P + g) 2 ' 

where the definition of B s meson decay constant is 

2 

(B s \bi l5 s\0) = (24) 

Combining PU, ([MD and (|23J), we have 

(mi) + m s )[m^ — (P + q) z \ °o ^ ^ s o 
°° , p h + (s,q 2 )P^ + p h _(s,q 2 ) q ^ 

S "(^ + g) 2 ' (25) 

where we have expressed the contributions from higher states of the B s channel in the form of 
dispersion integral with Sq b being the threshold parameter corresponding to the B s channel. 

On the theoretical side, the correlation function (|22|) can be also calculated in the perturbative 
theory with the help of the OPE technique at the deep Euclidean region P 2 ,q 2 = —Q 2 <C 0: 

U,(P,q) = n^ CD (q 2 ,(P + q) 2 )P,+U QCD (q 2 ,(P + q) 2 )q, (26) 

l ImnQ CD (q 2 ,(P + q) 2 ) „ , f 1 Imng CD (q 2 , (P + q) 2 ) 

(m b +m 3 ) 2 K S-(P + qy J{m b +m s ) 2 ^ S - {P + q) z 

Making use of the quark-hadron duality 

pHs,q 2 ) = hmIl? CD (q 2 ,(P + q) 2 )G(s-s h ), (27) 

with i = "+, — " and performing Borel transformation on both sides of Eq. (|27p with respect to 
(P + q) 2 , the sum rules for the form factors can be written as 

(m b +m s ) 2 "~ ^ ' Uy ^ V M 2 

To the leading order of a s , the correlation function can be calculated by contracting the bottom 
quark fields in Eq. (|22p and inserting the free b quark propagator 

U,(P,q) =ij ( d 4 xf ^2 ."fez (DU(PMxhMt + m b )i l5 s(0)\0). (29) 



2 m b + m s ttQCD/ 2 \ ( m B a ~ s x /nn\ 

/*(? ) = —; 2~ / dslmn^ (g ,s)exp( — 775— )• (28) 
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It should be pointed out that the full quark propagator also receives corrections from the back- 
ground field [3, Q] , which can be written as 



+ 2 1 12 vx^ u {vx) lu ], (30) 

where the first term is the free-quark propagator and G^J = G^ v T?j with Tr[T a T b ] = ^S ab . Sub- 
stituting the second term proportional to the gluon field strength into the correlation function can 
result in the distribution amplitudes corresponding to the higher Fock states of .0^(2317) meson. 
It is expected that such corrections associating with the LCDAs of higher Fock states do not play 
any significant roles in the sum rules for transition form factors 
neglected. 

Substituting Eq. (|17p into Eq. fl29f) and performing the integral in the coordinate space, the 
correlation function in the momentum representation at the quark level can be written as 

U,(P,q) = /'^[(ms+timDj^t^ (31) 

with 



351 ]. and hence can be safely 



(P 2 + q 2 - s ) + J(P 2 +q 2 - s ) 2 + 4P 2 K " Q 2 ) 

^ = " W 2 • (32) 

Combining Eq. (|28|) and Eq. (|3ip . we can finally derive the sum rules for form factors (q 2 ) 

s0 m B a JB s Juo u 

B. Sum rules for B® — > £)+(1968) transition form factors 

The form factors responsible for the B® — > Df transition are defined by 

(Df(P)\c^b\B s (P + q)) = fi s (q 2 )P, + fB s (Q 2 )^- ( 34 ) 
The correlation function associated with the form factors fp (q 2 ) and fp B (q 2 ) can be chosen as 

MP,q) = ~ J d 4 x e ^(D+(P)|T{i 2M (x) ) j 1 (0)}|0), (35) 
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where the current 32a{%) is given by 



hjx = c(x)7 At fe(x). 



(36) 



One can write the phenomenological representation of the correlation function at the hadronic level 
simply by repeating the procedure given above as 



(m b + m s )[m 2 Bs - (P + 



■ifUm + fBSn] 



(37) 



8- (P + q) 2 

On the other hand, the correlation function at the quark level can be calculated in the framework 
of perturbative theory to the leading order of a s as 

" x dui [(m b + um Ds )(j) Ds (u)P^ + m Ds cp Ds («)<7 M ] 



n„ (P,q) 



(38) 



iuo u s- (P + q) 2 

where uq has been defined in Eq. (|32|) . Matching the correlation function obtained in the two 
different representations and performing the Borel transformation with respect to the variable 
(P + q) 2 , the sum rules for the form factor f Bg (q 2 ) and f^ a (q 2 ) can be derived as 



tiM 2 ) 



m b + m s m 2 B JM 2 f 1 ^ 

m ljB s e J u 



— (m 6 + nm D J^(n)e-K-^ 2 +^ 2 )/M^), 
u 



{q2) = m + r^ emk/M , 

™>bJb. J u 



du 
u 



m Ds 4> D ( U )e-^-^ 2+uUp2 y^ M2 \ 



(39) 



IV. NUMERICAL ANALYSIS OF SUM RULES FOR FORM FACTORS 



Now we are going to calculate the form 
parameters used in this paper [g, 32 



36 



37 



"actors fn s (q 2 ) and f l D (q 2 ) numerically. The input 
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39, 



40J are collected as 



\V US \ = 0.226, 



\V CS \ = 0.973, 
1.20±° o f 7 ) x 10" 3 , \V ts \ = 40.96 x 10~ 3 , 



V ub \ = (3.68ig:J|) x 10~ 3 , 
V cb \ = (43.991°$) x 10 ~ 3 
Vft\ = 0.999, 



\V ud \ = 0.974, 
\Vcd\ = 0.225, 
\Vtd\ 

a = (90.6l 3 J)°, /3= (21.581° jj)°, 

m b = (4.8 ± 0.1)GeV, m c (lGeV) = 1.275GeV, m s (lGeV) = 142MeV, 

m Bs = 5.368GeV, m Ds = 1.968GeV, m D * o = 2.318GeV, 

f Bs = (151 ± 12)MeV f Ds = (273 ± 10)MeV , f D * Q = (225 ± 25)MeV . 



7 = (67.8l*|)°. 



(40) 
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38j] and fsJfB 



As for the decay constant of B s meson, we use the results fs = 130MeV 
1.16±0.09 [39( determined from QCDSR. The leptonic decay constants of £> s (1968) and L> s0 (2317) 



are borrowed from Ref. 



32, 



391 ] . The threshold parameter sq can be determined by demanding the 
sum rule results to be relatively stable in allowed region for Borel mass M 2 , and its value should 
be around the mass square of the first excited states. As for the heavy-light systems, the standard 
value of the threshold in the X channel would be s x = (mx + Ax) 2 , where Ax is about 0.6 GeV 



41 



42 



43 



44j, 



45 



46[ |) and we simply take it as (0.6±0.1) GeV corresponding to s Bs = (36±2)GeV 



for the error estimation in the numerical analysis. 

With all the parameters listed above, we can proceed to compute the numerical values of the 
form factors. The form factors should not depend on the the Borel mass M 2 in a complete theory. 
However, as we truncate the OPE up to leading conformal spin for the distribution amplitudes 
of D s meson in the leading Fock configuration and keep the perturbative expansion in a s to the 
leading order, a manifest dependence of the form factors on the Borel parameter M 2 would emerge. 
Therefore, one should look for a working "window" , where the results only vary mildly with respect 
to the Borel mass, so that the truncation is acceptable. 

In the first place, we focus on the form factors at zero momentum transfer. As for the form 
factor /£) (0) associated with B s — > D s transition, we require that the contribution from the 
higher resonances and continuum states should be less than 30 % in the total sum rules and the 
value of /o (0) does not vary drastically within the selected region for the Borel mass. In view 
of these considerations, the Borel parameter M 2 should not be too large in order to insure that 
the contributions from the higher states are exponentially damped as can be observed form Eq. 
(|39p and the global quark-hadron duality is satisfactory. On the other hand, the Borel mass could 
not be too small for the validity of OPE near the light-cone for the correlation function in the 
deep Euclidean region, since the contributions of higher twist distribution amplitudes amount to 
the higher power of 1/M 2 to the perturbative part. In this way, we indeed find a Borel platform 
M 2 G [8, ll]GeV 2 as plotted in Fig.Q] The value of /^(O) is 0.86t°;:i£, where we have combined the 
uncertainties from the variation of Borel mass, the fluctuation of threshold value, the uncertainties 
of quark masses and the errors of decay constants for the involved mesons. Following the same 
procedure, we can further compute the other form factors numerically, whose results have been 
grouped in Tabled 

Now, we can investigate the q 2 dependence of the form factors f}j s (q 2 ) and f l D *^(q 2 ). It is known 
that the OPE for the correlation function (|22l [34"|) is valid only at small momentum transfer region 
< q 2 < (nib ~ itt-c) 2 — 2Aqcd( ?7 ^6 — "ic)- As for the case with the large momentum transfer (small 
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8 8.5 9 9.5 10 10.5 11 2 4 6 8 10 

M 2 (GeV 2 ) (a) q 2 (GeV 2 ) (b) 



FIG. 1: Left figure: The dependence form factor (0) on the Borcl mass M 2 (red solid line), the contribu- 
tion of higher states in the whole sum rules (magenta short dashing line) and the contributions of the twist-3 
light cone distribution amplitude's contribution(blue long dashing line); Right figure: The dependence of 
form factor (q 2 ) on the momentum transfer q 2 within the whole kinematical region. 



recoil region), it is expected that HQET works well for the b — ► c transition. In the framework of 



HQET, the matrix elements responsible for B s — ► D sx transition can be parameterized as 



( D *so{P)\ci^5b\B s (P + q)) = -i^mB a m D *Jrfo.Jw){v + u% + rj D *Jw)(v - v%], 

(Df(P)\c^b\B s (P + q)} = Vm B .mD,ftW(» + + r,^(w)(v - v%], (41) 

where v = (P + q)/m<B s and v' = P/m,D ax are the four- velocity vectors of B s and D sx mesons with 
D sx being either D s or D* Q meson and w = v ■ v', Combining Eqs. (J2Tj) . (|3~4"j) and (|4ip. we have 

1 



:[(m Bs + m Dsx )r]J{w) - (m Ba - m Dsx )Vi H], 



(42) 



y/m Bs m D 
V m B s 

with w = (m Bs + fn 2 £) ax — Q 2 )/^ m B s m D 3X - I n the heavy quark limit, the form factors rj^ (w) and 
Vd s ( w ) satisfy the following relations 



(43) 



where £(w) is the Isgur-Wise function 47|, |48| with the normalization £(1) = 1. Similarly, heavy 
quark symmetry allows to relate the form factors r/J* (w) and (w) to a universal function 
t 1/2 {w) [49| 



Vd* (w) + V D * ( W ) = - 2t i/2(w), Vd* b0 ( w ) - V D * s0 (w) = 2ti/ 2 (w). 



(44) 
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An important relation between the B — > D** form factors at zero recoil region and the slope p 2 of 
the B — > D^*> Isgur-Wise function is 



n m 

under the name of the Bjorken sum rule j^. Here, D** denotes the generic L = 1 charmed states, 



the subscript n, m identify the radial excitations of the states with the same J p . For the B -> D 
transition form factors, the essential difference with the Isgur-Wise function £(y) is that one can 
not invoke heavy quark symmetry arguments to predict the normalization of Ti^iw) 

Phenomenologically, one can parameterize the B s — ► D s (1968, 2317) form factors in the small 
recoil region as 



r]f(w) = r,f{l) + af(w - 1) + bf(w - l) 2 , (46) 

where the r)f(w) denotes the form factor rj D (w) and 77™ (w). The parameters r]f(l), af and 
bf can be determined under the condition that the form factors derived in the LCSR and HQET 
approaches should be connected in the vicinity of region with q 2 ~ (mj, — m c ) 2 — 2AQco( r ^b — Tn c ). 
In this way, we can derive the results of form factors in the whole kinematical region as shown in 
Fig. ([T]) as an example. The values of all form factors are tabulated in Table HI where the results 
under the QCDSR approaches are also collected for comparison. 

As can be observed from Table HI the number of form factor (w) at the zero-recoil region 
deviates from the zero significantly indicating that the l/m c power correction is sizeable for the 
B — > D* transition. Generally, the expansion of the current 

cTib = c V2 Tib Vl - -^—c V2 Tiip 2 b Vl + 7: — c^T^M^i + ■■• (47) 

constitutes the main source of the power corrections. The QCDSR estimation of the form factor 
/n* differs from that obtained in the LCSR approach in sign implying that the power corrections 
and radiative corrections of correlation function are in need to reconcile the existing discrepancy 
between these two methods. 



V. SEMILEPTONIC AND NONLEPTONIC DECAYS OF B s -> D s (1968, 2317) 

With the form factors derived above, we can further investigate the semileptonic and nonleptonic 
decays of B s to D s (1968, 2317) in the factorization approach. Factorization theorem is a basic 
theoretical tool to disentangle physical effects from different energy scales. Factorization in heavy 
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TABLE I: Numbers of /1(0) and r]f(w) determined from the LCSR approach, where the uncertainties from 
the Borel mass, threshold value, quark masses and decay constants are combined together. For comparison, 
the results estimated in the QCDSR are also collected here. 



this work QCDSR 



f+Jqi) 0.53t°;i? 0.40 ±0.10 [5] 

f D Jq 2 ) 0.18l°;°| -0.12 ±0.13 [5] 

/+( g 2 ) 0.86t°;£ 0.70[50] 

f Ds (q 2 ) 0.26±g:8| 0.19[50] 



^(1) 



„+» 0.45^;^-0.79l°;^ 0.92lg;g 

0.0518S -o.i2±S:8S o-iTlg^ 
H 0.7818:2 -i.i2l°i° o.86±g : i2 

VDs H 0.091^-0.161^ 0.14±8:S 



quark decays was firstly proposed in Ref. 5l[ as a phenomenological approximation similar to 
the vacuum saturation approximation for the four-quark operator matrix elements. Intuitively, 
factorization might work at least to leading-order approximation, since the partons that eventually 
form the emitted meson escape from the heavy meson remnant as an energetic, low mass, color 
singlet system, therefore independently from the remnant. 



A. Semileptonic decays of B° s -> D+(1968, 2317)ZP; 



With the free quark amplitude given above and the transition form factors derived in the LCSR, 
we arrive at the differential decay width for B® — > D+(1968, 2317)/P; modes 



dl 

dq 1 



G 2 F \V cb \ 2 (q 2 



m 



2\2 



7687T 3 



m 



B 



(q 2 ) 3 



, a 2 2/ 2 

+bq m l [m B 



m L>x 



-^(2m 2 (X + 3 q 2 m 2 Dx ) + q 2 X)\ft(q 2 )\ 2 
-Q 2 )ft(Q 2 )fr(Q 2 ) + ^rn 2 \fr( q 2 )\^ 



(48) 



m 



,2\2 



2™2 



Dx 



Aq m 



Dx- 



with A = {m 2 B 

For convenience, the q 2 dependence of these invariant functions are also plotted in Fig. [2] and 
[HJ Integrating Eq. (|48p . we get the branching fractions of — ► D+(1968, 23l7)Wi as grouped 
in Table [III It can be observed from this table that the orders of magnitudes for BR(B® — ► 
I?! (1968, 2317)Zi^) obtained in the quark model and sum rule approaches are consistent with each 
other. Besides, we can also find that the decay rates for the final state with r lepton are generally 
3 — 4 times smaller than those for the muon case due to the suppression of phase spaces. Once the 
data on the B® — ► Df (2317)Wi are available, the theoretical predictions presented here can be put 
to the experimental scrutiny to test the ordinary cs picture of -D s (2317) meson. 
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B. Nonleptonic decays of B s -> L>+(1968, 2317)M 

Now, we turn to the calculations of nonleptonic decays B s — > _D+(1968, 2317)M, where M can 
be a light meson or a charmed meson. As mentioned above, the factorization assumption will be 
employed to decompose the matrix element of four-quark operator 

(D SX M\Q\B S ) = (Mlh^iDsMBs), (49) 

into the the B s — ► (1968, 2317) transition form factors and the decay constant of M. 

For a light meson M, only the tree-operators in Eq. ([2]) can contribute to these decay modes 
induced by the b — > c transition. Then, the decay width for B s — > L>^"(1968, 2317)L can be written 
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TABLE II: Branching ratios for the semilcptonic decays — > D+(1968, 2317)Zi^ with the form factors 
estimated in LCSR, where the results calculated in constituent quark model and QCDSR are also displayed 
for comparison. 

B°^D+ 1-Pi l = e,fi 1 = t 

this work (2.3tJ;o) x 10 " 3 ( 5 -7±||) x 10~ 4 

QCDSR[5J ~ 10" 3 ~ 10~ 4 

Constituent Quark Model[3j (4.90 - 5.71) x 1CT 3 
QCDSR in HQET[4] (0.9 - 2.0) x 10~ 3 

B° s ^D+l-Di l = e,fi 1 = t 

this work (l-Otol) x 10~ 2 (3.3±l;f) x 10" 3 

Constituent Quark Model[3] (2.73 - 3.00) x 10~ 2 
QCDSR [6] (2.8-3.8) x 10- 2 



as 



Gl\P\ n f \fn (mj)m L (e* ■ P)\ 2 (L=V), 



16vr< lc - l r B -->- m V ^K) + /^(m|)m|P (L=5, 

where L denotes a light meson; V, P and S* label the vector, pseudoscalar and scalar mesons 
respectively. The magnitude of the three-momentum for the recoiled charmed meson is 
,p, = i( m B s ~ ( m Ds X + m L ) 2 ){m 2 Ba - {m Dsx - m L ) 2 )} 1 / 2 

' ' 2m Bs ' 1 ' 

and the decay constant Jl has been collected in Table lllli The decay constants for the light 
pseudoscalar mesons are taken from the Particle Data Group 

decay constants are extracted from the data on r~ — ► (p~ ,K*~)i/ T . To determine the decay 
constants for the scalar meson Dq and vector meson D*, the following relation 



361 ] and the vector meson longitudinal 



Id* Id* t 



B 



(52) 



is assumed in this work. The decay constant of vector meson D* is borrowed from Ref. [39(. The 
energy scale of the Wilson coefficient ci2(/-0 is varied from 0.5m;, to 1.5m;, in the error estimations. 

Substituting the form factors obtained in the previous sections into Eq. (I50j) . we can get the 
decay rates of B s — > Z?+(1968, 2317)L as shown in Table ITVl From this table, the results evaluated 
in the factorization approach are in accord with that predicted in the PQCD approach and the 
available data, which implies that the factorization assumption (FA) works well for these color 
allowed modes as expected. 
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TABLE III: Decay constants of light and charmed mesons (unit: McV). 



u 


Sk 


fp 


!k> 


So 


fo- 


f D - 




131 


160 


209 ±2 


217±5 


206 ±9 


95 ± 10 


270 ± 35 


312 ±40 



TABLE IV: Branching ratios(unit: 10 -4 ) for the nonleptonic decays B s — > Df x L (L denotes a light meson) 
estimated in the FA with the form factors obtained in the LCSR, where we have combined the uncertainties 
from the form factors, the scale-dependence and CKM matrix elements. In Rcf. [50j . the authors also 
employ the naive factorization but take the transition form factors calculated in the three-point sum rules. 



Channels 


this work 


PQCD[JJ 


Exp. [36] 


FA [50] 


B° s ~ 












B° s - 


* D^K- 










B° s ~ 


- D%p~ 


13+ 6 5 








B° s ~ 


> D*+K*- 










B° s - 


-> D+TT- 


i '-6 


1Q c + 10.6+6.3+0.6 
l».U_ 7 5 _g .2-0.6 


32 ±9 


11.2 ±22.7 


m- 


-> D+K- 


i o+0.5 
i -°-0.4 


i 7n +0. 87+0. 53+0. 05 
i - ' u -0. 66-0. 56-0. 05 




0.83 ±2.51 


B° s 


- D+p- 


42 +i7 


A7 n +24.9+15.3+1.36 
+ ' ^-17.7-14.8-1.37 






m- 


+ DfK*- 


2 4 +1 '° 


9 81 +1-47+0.79+0.09 
z - 01 -l. 09-0. 85-0. 09 




0.60 ±2.82 



Moreover, it is also helpful to define the following ratios 
BR{m - D%ir~) BR(B°^D+7r- 



Ri 
R2 



BR{B>^D*+p-) ^ BR(B° s ^D+p-) 



Vud fn 
Vus fl< 

Vud fp 



V us f 



A * 



12.6, 
« 17.4, 



(53) 



5i2(.B° D*q A*~) BR(B° DJK*-) 

which are consistent with those collected in Table HVl 

As for the two charmed meson decays of B s meson, the decay width in the factorization approach 
can be given by 



F(B S -> D+X) 



G%\P\ 



\V cb V* q a 2 {p) - V tb V tq [a A {p) + a w (p) + r q (a 6 (p) + a 8 (/x))]| /. 



2 t 2 
X 



|/+ (ml)m x (e* ■ P)\< 



-f£j m x) + f D J m x) m 



x 



(X=V), 
(X=S, P), 



(54) 
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with 



2m 2 x 




-m c )(m c +m q ) 




2m 2 x 


K 


—m c )(m c —m q ) 




2m 2 x 


(m b - 


m c )(m c -m q ) 




2m 2 x 


K- 


m c )(m c -in q ) 












=P), 




- X-- 


=S), 




> D%X, X-- 


=P), 




* D*+X, X-- 


=S), 


(5?- 


- D+X, X-- 





(55) 



where the quark masses in the above equation are the current quark masses. 

Combining the Eq. (]54j) and the form factors listed above, one can easily get the branching 
ratios of B s — > D+.X (X being a charmed meson) as shown in Table [V] As can be seen from 
this table, the decay model B® — ► D+ D~ possesses a quite large branching ratio of order 10~ 2 , 
which should be detected easily at the large colliders such as Tevatron and LHC. Moreover, the 
theoretical predictions on the B® — ► D+tt~ decay can accommodate the experimental data within 
the error bars. As for the B® — > D+D*~ decay, only the upper bound for this mode is available at 
present, which is also respected by our predictions. 

Subsequently, the ratio of decay rates between the Cabibbo favored and suppressed modes can 
be estimated as 

BR(B® 



R 3 = 
-R4 = 
R5 = 



BR{B 



D*+ D~ 



D+D: 



D*+D~) BR(B° - D+D-) 



V cs fn- 



BR{B 



D*+ D* 



BR(B° S 



DfD* s 



Vcd /d- 



BR(Bg 



D*+D*~ 



BR(B* 
BR(B° S 



D+D 
D+D 



sO 



r)*+ f)*- 



D+D 







Vcd Id*- 

V cs fD*n 



v cd f 



IK, 



(56) 



BR(B$ 

in the naive factorization without the contributions from penguin operators. Such naive estimations 
are in good agreement with that presented in Table |Vj which also indicates that the two charmed- 
meson decays of B s meson governed by the b — > qcc (q = s d) transition are dominated by the tree 
operators. 



VI. DISCUSSION AND CONCLUSION 



A detailed analysis of properties about the new charming mesons such as D s (2317) has become 
a prominent part of the ongoing and forthcoming experimental programs at various facilities world- 
wide. The production characters of charmed mesons in the B s decays are especially interesting for 
highlighting the understanding of QCD dynamics and enriching our knowledge of flavor physics. 
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TABLE V: Branching ratios(unit: 10 -3 ) for the nonleptonic decays B s — > Df x X (X denotes a charmed 
meson) estimated in the FA with the form factors obtained in the LCSR, where the uncertainties from the 
form factors, the scale-dependence and CKM matrix elements have been combined together. 



Channels This work Exp. [36] 







I3±l 




m- 




o.5l8:2 




m- 


- DSD* 


^•■■—0.8 




b° s - 








m- 




6.oll:S 




m- 


> D*+D*~ 


0.2l° j 




b° s - 


- DfD- 


35^ 


110 ±40 


B° s - 


- D+D- 


1 l +0 - 4 

1 - 1 -0.4 








c o+2.2 
J - J -1.8 




m- 


- DfD*- 


o. 2 ±° :l 




m- 




331}? 


< 121 


B° s - 




i-4l8:S 





More importantly, the theory underlying the description of the decays induced by the b — > c tran- 
sition is mature currently. In view of the large mass of b quark, heavy quark expansion works well 
enough to enable a precise determination of the decay amplitude. 

LCSR approach is employed to compute the B® — > Df (1968, 2317) transition form factors at 
the large recoil region and the results are then extended to the small recoil region in the framework 
of HQET. Our results show that the power correction to the form factor r^, (w) responsible for the 
£° -> £>+(1968) at the zero-recoil region transition is numerically small, since this form factor only 



52J. However, 



receive the corrections at order 1/m^ and l/m^ as indicated by the Luke's theorem 
the power correction to the form factor (w) relevant to the B — > D*^ transition is sizeable. 

Subsequently, we utilize the form factors estimated in the LCSR approach to perform a careful 
study on the semileptonic decays B® — ► -D+(1968, 2317)/P/. It has been shown in this work that 
the branching fraction of the semileptonic B® — ► D+(2317)^P^ decay is around 2.3 x 10~ 3 , which 
should be detectable with ease at the Tevatron and LHC. The decay rates of semileptonic modes 
for the final states with r lepton are approximately 3 — 4 times smaller than those with muon due 
to the suppression of phase spaces. In addition, the branching fractions of B® — ► L>+(1968)Zz^ are 
almost one order large than that for the B® — > Df (23l7)Wi decays. 
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Nonleptonic decays B s — > Df(1968, 2317)M are also investigated in the framework of factoriza- 
tion approach in this work. It is found that the theoretical predictions for B s — > D+(1968, 2317)L 
presented here are in agreement with those obtained in the kx factorization, which supports the suc- 
cess of color transparence mechanism in the color allowed decay modes. Moreover, B® — > D*^Dj 
owns a quite large branching ratio as 1.3 x 10~ 2 , which should be accessible experimentally. More 
theoretical results worked out here are expected to be tested by the large colliders in the near 
future. 
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